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ABSTRACT 

The pressure profiles of the Intracluster Plasma in galaxy clusters show a wide variance when 
observed in X rays at low redshifts z < 0.2. We find the profiles to follow two main patterns, featuring 
either a steep or a shallow shape throughout both core and outskirts. We trace these shapes back 
to a physical dichotomy of clusters into two classes, marked by either low entropy (LE) or high 
entropy (HE) throughout. From X-ray observations and Sunyaev-Zel'dovich stacked data at higher 
0-2 ^ z ^ 0.4, we elicit evidence of an increasing abundance of HEs relative to LEs. We propose 
this to constitute a systematic trend toward high z; specifically, we predict the pressure profiles to 
converge into a truly universal HE-like template for z > 0.5. We submit our physical templates and 
converging trend for further observational tests, in view of the current and upcoming measurements 
of individual, stacked, and integrated Sunyaev-Zel'dovich signals. 

Subject headings: cosmic background radiation — galaxies: clusters: general — X-rays: galaxies: 
clusters — methods: analytical 



1. INTRODUCTION 

A keen interest is focusing on the radial pressure pro- 
files p{r) that prevail in the hot intracluster plasma (ICP) 
filling up the galaxy clusters. Specifically, it is debated 
the degree of their 'universality' among rich clusters, on 
the following grounds. 

On the upside, the profile of the ICP thermal pres- 
sure p = nkBT/fj, (with the mean molecular weight 
/i « 0.59 for cosmic abundances) is expected to score off 
the components: temperature T(r) and number density 
n{r) = p[r)/mp, on account of its prompt equilibration 
at sound speed in the absence of forcing stresses. The 
equilibrium gradient is simply linked by 



dp 
d7 



GM{< r) 



(1) 



to the gravitational force from the dark matter (DM) 
distribution M(< r) on the ICP mass density p. 

On the downside, one may worry that the equilibrium 
pressure profiles might be affected by the complexities 
elicited for redshifts z < 0.2 in the X-ray observations. 
These probe n from the surface brightness Sx oc n? T^/^ 
emitted by the ICP via thermal bremsstrahlung, while 
the temperatures fcsT in the keV range are measured 
from spectroscopy. 

In fact, in the cluster cores for r < 0.2i?5oo ^ the ra- 
diation observed in X rays is to erode the ICP thermal 
content on the cooling timescale tc ~ 30 (fcsT/keV)^/^ 
(n/lQ-^ cm-3)-i Gyr (e.g.. White & Rees 1978; Voit 
& Bryan 2001). The process tends to speed up as n 

^ We adopt the flat cosmology with matter density parameter 
= 0.3, Hubble constant Hq = 70 km s~^ Mpc~^, and mass 
variance erg = 0.8 (Komatsu et al. 2011). The virial radius reads 
-R ~ -Rioo ~ 4iJ20o/3 ~ 2 H500 in terms of the radii encircling 
average overdensities 100, 200 and 500 over the critical density, 
and takes on values around 2 Mpc for rich clusters. 



raises and T lowers, to the effect of steepening the inner 
pressure gradient after Eq. (1). A catastrophic runaway 
is conceivably offset by the energy fed back from self- 
regulated AGN activities (see discussions by Cavaliere et 
al. 2002; Lapi et al. 2003, 2005; Ciotti & Ostriker 2007; 
Churazov 2010), or is even forestalled by strong energy 
injections from deep mergers (see McCarthy et al. 2007; 
Markevitch & Vikhlinin 2007). 

In the cluster outskirts for r > R500, a pressure jump 
is forced by shocks driven by the supersonic inflow of 
cold external gas. Thus the gravitational infall energy 
is mainly thermalized at the virial boundary (see Lapi 
et al. 2005, 2010; Voit et al. 2005); meanwhile, part of 
it drives outer subsonic turbulence contributing to the 
equilibrium so as to require lower thermal pressures (see 
Lau et al. 2009; Cavaliere et al. 2011a). 

Bypassing such complexities, a 'universal' fitting for- 
mula for the pressure profiles has been proposed by Nagai 
et al. (2007) to interpret the outcomes of hydrodynami- 
cal simulations of relaxed clusters (see also Battaglia et 
al. 2011), and applied by Arnaud et al. (2010) to ren- 
der with empirically adjusted parameters the X-ray data 
out to i?5oo for z < 0.2. Actually, these analyses led 
to recognize an average profile along with a considerable 
variance. 

On the other hand, the pressure profiles can be di- 
rectly probed with the thermal Sunyaev-Zel'dovich ef- 
fect (1980; SZ) that occurs as CMB photons are in- 
verse Compton scattered by the hot ICP electrons, and 
change the radiation temperature Tcmb ~ 2.73 K by 
an amount AT = guyTc-mh ~ —0.5 mK. This provides 
a linear, intrinsically z-independent probe of the ther- 
mal electron pressure Pe = Mp/^e ~ 0.52 p (with the 
mean molecular weight per electron /ie « 1-15), since 
its strength is given by the Comptonization parameter 
y = {uT/meC^) J d£pe{r) ^ 10~* integrated along the 
1.0. s. The spectral factor g^, approaches the value —2 at 
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low frequencies (see Rephaeli 1995); its positive signa- 
ture for V > 217 GHz offers a powerful cross-check for 
the SZ nature of the signals. 

The SZ observations now start to probe the radial pro- 
files in nearby individual clusters, and in more distant 
stacked samples {SPT collaboration, Plagge et al. 2010; 
WMAP collaboration, Komatsu et al. 2011; Planck col- 
laboration, Aghanim et al. 2011a, 2011b). They arc also 
addressing the cluster contribution to the CMB power 
spectrum at multipoles i > 2000 (see Lueker et al. 2010; 
Dunkley et al. 2011; Reichardt et al. 2011). Extensive 
data at higher resolutions and sensitivities are expected 
from current and upcoming instrumentation, and even- 
tually from ALMA (sec Birkinshaw & Lancaster 2007). 
All such actively pursued observations call for a reliable 
template to interpret and assess their astrophysical and 
cosmological import. 

Here we take advantage of the effective formal- 
ism provided by the Supermodel (SM; Cavaliere 
et al. 2009, Fusco-Femiano et al. 2009, see 
http : //people . sissa . it/~lapi/ Supermodel/) 
to show that the complex ICP thermal states still allow 
a neat physical description of the spherically-averaged 
pressure profiles. In response to the intriguing challenge 
posed by Arnaud et al. (2010), we find two basic shapes 
to span the observed variance in the pressure profiles for 
low z < 0.2, and predict that they are to converge into 
a closely universal one for high z > 0.5. 



2. A PHYSICAL APPROACH TO PRESSURE PROFILES 

In singling out the ICP disposition and evolution we 
base on the updated paradigm for the hierarchical forma- 
tion of the containing DM halos (e.g., Zhao et al. 2003; 
Genel et al. 2010; Wang et al. 2011). The paradigm 
comprises an early collapse punctuated by major merg- 
ers, building up the core over a few crossing times; this 
occurs at redshifts Zf w 1.5 — 0.5 weakly depending on 
the mass M ~ 10^'' - 10^^ Mq. Slow, dwindling inflows 
of external matter follow over several Gyrs, and accrue 
the outskirts out to the current virial R in closely sta- 
tionary conditions described by the Jeans equation. In 
the process, the core scale radius r_2 stays put while R 
expands, and the 'concentration' parameter c = R/r-2 
correspondingly grows. 

The ICP forms from intergalactic gas with pressure < 
IQ-^ eV cm-3 (see Ryu et al. 2008; Nicastro et al. 2010) 
that - along with the DM - inflows at supersonic Mach 
numbers A4 > 10. The gas is shock-heated at about R, 
and its pressure jumps by factors 4 A^^ > 500 up to the 
ICP values pn > 1 cV cm^'\ Inward of R, the pressure 
rises to balance the DM gravitational pull as described 
by Eq. (1), with the following universal features: the rise 
will be monotonic within a smooth potential well; the 
gradient will vanish at the center when the gravitational 
force does, as implied by the Jeans equilibrium and found 
in simulations and real data (see Lapi & Cavaliere 2009a, 
2009b; Navarro et al. 2010; Newman et al. 2011); the 
rise terminates with a flnite central value, proportional 
to the thermal energy density. 

Within such universal constraints, a full description of 
the pressure profile is keyed to the specific 'entropy' (adi- 
abat) run k{r) = ksT/n'^^^ embodying the ICP thermal 



state. In fact, using n oc (p/k)^^^ Eq. (1) solves to yield 



p{r) 



mpGM{< x) 

J.2 /(;3/5^2;) 



5/2 



(2) 



in the context of the SM (see Lapi et al. 2005; Cavaliere 
et al. 2009). For the DM mass distribution M{< r) 

we use our a— profiles; these solve the Jeans equations 
under physical boundary conditions, and agree well with 
A''— body simulations (see Lapi & Cavaliere 2009a). 

The key role is played by the entropy run fc(r) over 
scales larger than some 10^ kpc, for which physical in- 
sight is provided by the above picture of cluster forma- 
tion. The latter indicates the basic shape (see also Voit 
2005) 

k{r) = ke + kH {r/RT ; (3) 

this involves the boundary scales R and fc/j, and features 
two intrinsic parameters: the central level fee, and the 
average outer slope a. These arc evaluated as follows. 

In the cluster core r < 0.2 i?5oo a level kc ~ 10^ keV 
cm^ is set by the early collapse; this results from den- 
sities n ^ 10~^ cm~^ compressed by the standard con- 
trast factor 200 over the average background's, and from 
temperatures impulsively raised to values /csT ^ GM(< 
r)/10r ^ a few keVs. Thereafter, radiative cooling com- 
petes with energy injections from AGNs or mergers to the 
effect of stabilizing or even raising the time- integrated kc 
at levels that gather around 10^ or 10^ keV cm^ (see Cav- 
agnoloetal. 2009; Pratt et al. 2010; Hudson et al. 2010). 
From Eqs. (1) and (2) the central pressure is seen to fol- 
low the basic scalings Pc oc kc^^^ and dp^/^/dr a kc^^^. 

At the other end r > i?5oo, an entropy ramp rises with 
slope a < 1.1, originated by the continuously shocked 
infall and the progressive stratification of the accreted 
shells during the slow outskirts growth (see Tozzi & Nor- 
man 2001; Lapi et al. 2005). The ramp ends up at 
the bounday R with the value fc^ > 10^ keV cm^ set by 
strong shocks (see Cavaliere et al. 2011b). Prom Eqs. (1) 
and (2) the outer pressure profile follows the basic scal- 
ing p{r) ~ PR (r/i?)~^+^", see Cavaliere et al. (2011b). 
We stress that p{r) will decrease monotonically outwards 
even when the temperature T(r) oc p'^^^{r) k^^^{r) fea- 
tures a middle peak at r « 0.2i?5oo due to the entropy 
rising steeply from a low kc] the peak is the defining mark 
of cool-core clusters (see Molendi & Pizzolato 2001; Lec- 
cardi et al. 2010). 

When the X-ray data on brightness and temperature 
of clusters at z < 0.2 are analyzed with the SM, a direct 
correlation between kc and a emerges (sec Cavaliere et 
al. 2011b). This implies that clusters can be parted 
into two main classes: LE or HE, featuring low or high 
entropies, respectively, throughout cores and outskirts. 
The LEs (e.g., A2204, A1795) are marked by low central 
entropies kc ^ 10^ keV cm^ along with shallow entropy 
slopes a < 0.7; the HEs (e.g., A1656, A399) are marked 
by higher central values kc > 10^ keV cm^ along with 
steeper slopes a ~ 1. 

Then from the above inner and outer scalings 
dp^/^/dr oc kc^^^ and p oc t-"^"""^", we find the pres- 
sure profiles to differ from HEs to LEs, with the former 
featuring quite shallower gradients both in the core and 
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Figure 1. Profiles of ICP pressure normalized to psoOi see § 3 for details. The yellow shaded area illustrates the region covered by the low 
redshift (z < 0.2) clusters of the REXCESS X-ray sample; the dotted blue and red lines refer to the average profiles for the subsamples of 
cool-core and non-cool-core clusters, separately. The cyan shaded area illustrates the coverage by hydrodynamical simulations of relaxed 
clusters. The dashed line represents the joint fit to the observational and virtual data with the 'universal' pressure profile. The green 
squares represent stacked SZ observations of higher redshift (0.2 < 2 < 0.4) clusters with the SPT. Our SM templates for HE and LE 
clusters are illustrated by the red and blue solid lines, respectively. 



in the outskirts. Our picture is substantiated by the data 
presented in Fig. 1. 

3. PROBING PRESSURE WITH X RAYS AND SZ EFFECT 

In Fig. 1 we compare the pressure profiles computed 
with our SM to data from X-ray and SZ observations, 
and from numerical simulations. In the plot, the ra- 
dial scale is normalized to -R500, while the pressure 
is normalized to the standard value p^oo ~ 1.8 /iz^^ 
(M500/5 X 10"Mq)2/3 eV cm-3 in terms of the Hub- 
ble parameter = H{z)/Hq (e.g., Ettori et al. 2004; 
Arnaud et al. 2010). 

Our pressure templates provided by Eq. (2) for HE 
and LE clusters are illustrated by the red and blue solid 
lines; these are computed from Eqs. (2) and (3) with the 
parameter values discussed in § 2. Specifically, for typical 
HEs we adopt = 100 keV cm^, fc^. = 3 x io^ keV cm^ 
and a = 1.1, while for LEs we adopt kc — 10 keV cm^, 
kji — 10^ keV cm^, and a = 0.7; for both we set i? = 2 
Mpc. These values are consistent with the outcomes from 
detailed SM fits to X-ray observations of nearby clusters 
carried out by Fusco-Femiano et al. (2009) and Cavaliere 
et al. (2011b); note, however, that fitting the centrally 
flat profiles of HEs (e.g., A2256) requires the level kc to 
extend out to r/ ^ lO'^ kpc, bearing the imprint of a 
recent merger. 

The yellow shaded area illustrates the region covered 



by the low redshift (z < 0.2) clusters of the REXCESS 
X-ray sample analyzed by Arnaud et al. (2010); the dot- 
ted blue and red lines refer to their average profiles for 
the subsamples of cool-core and non-cool-core clusters, 
respectively. The cyan shaded area illustrates the region 
covered by hydrodynamical simulations of relaxed clus- 
ters (Borgani et al. 2004; Nagai et al. 2007; Piffaretti & 
Valdarnini 2008; Battaglia et al. 2011). 

The dashed line represents the joint fit by Arnaud et 
al. (2010) to the observational and virtual data in terms 
of their 'universal' pressure profile. For r < R5Q0 the 
X-ray data show that such a profile yields only an aver- 
age but incomplete description. In fact, the partial av- 
erages over the cool-core and non-cool-core subsamples 
deviate upward and downward by a large amount exceed- 
ing their internal variance; thus a bimodal description 
constitutes both a closer and a more effective represen- 
tation. This is just what is provided by the above SM 
templates for HE or LE clusters, which in the core recover 
the non-cool-core or cool-core behaviors. Moreover, for 
r > i?5oo where only scarce X-ray data are available, the 
SM template for LE clusters agrees well with the results 
of hydro— simulations of relaxed clusters. In the way of a 
prediction, beyond i?5oo we expect for HEs considerably 
higher pressure profiles relative to LEs, as represented in 
Fig. 1. _ 

In addition, our picture envisages for decreasing z an 
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Figure 2. Left panel: Profile of the SZ signal from the Coma Cluster, in terms of Rayleigh-Jeans temperature decrement AT. The green 
filled dots illustrate the WMAP data. The black dashed line refers the 'universal' pressure profile. The red and blue solid lines illustrate our 
SM outcomes for HE and LE clusters. Right panel: Cosmological thermal SZ effect at 150 GHz. The black dot with error bars illustrates 
the current constraints set with the SPT. The dashed line refers to the 'universal' pressure profile. The outcomes from our SM templates 
are illustrated by the solid lines: blue for LEs, red for HEs, and orange for HEs with a nonthermal component increasing with z (see § 3). 



ICP evolution from HE to LE states. In fact, over life- 
times of several Gyrs elapsed from the formation w 1 
to the observation redshift z < 0.2 the outer slope a will 
flatten out from values a « 1.1 towards values a < 0.5, 
and correspondingly the boundary level kfj decreases; 
such a trend is enhanced as z approaches 0, when the 
cosmic timescale lengthens considerably. This is because 
the entropy production by weakened virial shocks is re- 
duced as inflows peter out, especially in the accelerated 
cosmology (see Lapi et al. 2010). Meanwhile, cool- 
ing erodes the central entropy over comparable times 
tc « 5(fcc/100keVcm2)i-2 Gyr after Eq. (2), and un- 
dergoes an accelerated drop toward the 'attractor' level 
kc ~ 10 keV cm^ set by competition with AGN feedback. 

On the other hand, the evolution of the cluster DM 
halo is marked by the growth of the 'concentration' c w 
3.5 hzf/hz from Zf to z (see Zhao et al. 2003; Prada et al. 
2011). Since a and kc decrease together, we expect them 
both to anticorrelate with c. In fact, such correlations 
a, kc vs. have been quantitatively elicited with SM 
analyses of high-quality X-ray data concerning several 
clusters at z < 0.2 (see Gavaliere et al. 2011b). 

As z increases, our evolutionary trend envisages the 
HE/LE ratio to grow towards 1; in other words, we ex- 
pect all pressure profiles to converge toward a truly uni- 
versal template provided by the HE shape. Our pic- 
ture is supported by comparison of the local data with 
stacked SZ observations of redshifts 0.2 ^ z < 0.4 clus- 
ters; the pressure profiles from the SPT stacked data 
(Plagge et al. 2010) are represented in Fig. 1 with the 
green squares. Although the uncertainties are still con- 
siderable in the outskirts, a departure from the 'univer- 
sal' profile stands out, and the trend toward an HE- like 
template clearly emerges. The same trend is emerging 
from the analysis of a stacked cluster sample observed 
with WMAP (Komatsu et al. 2011). A similar trend is 
suggested by the sample of clusters detected by Planck 
for redshift 0.3 < z < 0.5, and followed up in X rays with 
XMM-Newton (Aghanim et al. 2011b). Independent evi- 
dence is provided by the dearth of strong cool-cores found 
in X rays by Santos et al. (2010) at high-redshifts. 



The next step to test the evolutionary trend from LEs 
to HEs will involve observing the SZ profile from in- 
dividual clusters over an extended range of redshifts. 
This is still challenging with the present instrumenta- 
tion even at z < 0.2, but will become feasible out to high 
z > 0.5 with new-generation instruments up to ALMA 
(see http : //www. almaobservatory . org/). To illustrate 
the current status, in Fig. 2 we report the recent data 
on the SZ profile for the Coma Cluster at z w 0.023 
with WMAP (Komatsu et al. 2011). We also report the 
'universal' pressure profile (black dashed line), and the 
Supermodel templates for HE and LE clusters (red and 
blue solid lines). These data allow to recognize the HE 
nature of the Coma cluster, still more prominent in X 
rays. On the other hand, the SZ data will become com- 
petitive once resolutions better than 1' will be attained. 

Another testbed for our evolutionary picture will be 
provided by the power spectrum of the unresolved SZ 
effect integrated over redshift and over the evolving clus- 
ter mass distribution including groups (e.g., Shaw et al. 
2010; Efstathiou & Migliaccio 2011). In Fig. 3 we il- 
lustrate the outcome from the 'universal' profile of Ar- 
naud et al. (2010) and the SM templates for LE and 
HE clusters, compared with the constraints at £ ^ 3000 
set by current observations with the SPT (Reichardt et 
al. 2011). The latter are converging to indicate that the 
integrated SZ effect is to be both dominated by HEs at 
the relevant z > 0.5, and reduced by the presence of a 
nonthermal contribution to the inner ICP equilibrium. 

We find consistency of our HE template with the data 

when the standard scaling psoo oc hl^^ is decreased by a 

— 1/2 

factor hz ; the latter renders how the inner nonther- 
mal component mildly increases with z, as expected for 
the early HE equilibrium punctuated by major mergers 
with strong turbulent wakes (see § 2 and Gavaliere et al. 
2011a). Additional constraints on the detailed shape of 
the SZ effect will require high sensitivities at resolutions 
£ > 5000 with good control of the systematics, an effort 
actively pursued by the SZ community. 
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4. DISCUSSION AND CONCLUSIONS 

As the observational uncertainties in the SZ measure- 
ments are shrinking below the theoretical ones, we have 
aimed at eliciting truly universal features in the pressure 
profiles of the Intracluster Plasma (ICP). These are de- 
scribed in terms of our Supermodel, are tested at low 
z against X-ray observations of brightness and tempera- 
ture, and are proposed for probing at higher z with the 
linear, z-independent SZ effect. 

At low z < 0.2 we identify from X-ray data two main 
pressure patterns, each featuring related shapes at the 
center and in the outskirts. These are interpreted in 
terms of two cluster classes: HEs with high entropy 
throughout, i.e., high central levels kc along with a steep 
outer rise with slope a, leading to shallow pressure pro- 
files; conversely for LEs. 

For z > 0.2 we expect the onset of an evolutionary 
trend in cosmic time from HE to LE, comprising fast 
erosion of central entropy by radiative cooling, along 
with reduced outer entropy production by decreasing in- 
flows (see Cavaliere et al. 2011b). Accordingly, the ratio 
HE/LE is to increase toward high z, implying for z > 0.5 
convergence of the pressure profiles toward a truly uni- 
versal, HE-like template. 

Few intermediate instances are expected; these may oc- 
cur when cooling is forestalled by an exceedingly high ini- 
tial level kc > 10^ keV cm^, while the outer entropy ramp 
is independently flattening (e.g., in A1689 and A2218 
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as discussed by Cavaliere et al. 2011b). On the other 
hand, Eq. (2) implies our radial pressure profiles to be 
robust against overall asphericity and localized hot/cold 
imprints from recent mergers, even when the latter cause 
wiggles in temperature and marked flatness in the cen- 
tral brightness as often found from detailed fits to X-ray 
data in HEs (e.g., A2256 as discussed by Fusco-Femiano 
et al. 2009). This strengthens the case for universality 
of our asymptotic HE-like pressure template. 

We have traced evidence for our evolutionary trend 
developing in the pressure profiles through stacked SZ 
signals due to 0.2 ^ z < 0.4 clusters (see Fig. 1). At 
higher z, further evidence will be difficult to pinpoint 
from X rays alone, given their bias toward high central 
brightnesses proper to LEs; rather, it will be provided 
again by SZ signals from individual or stacked clusters. 
In parallel, such an evidence will be tested with the in- 
tegrated contribution from unresolved SZ signals to the 
CMB anisotropics at muhipoles i > 5000 (see Fig. 2). All 
such observations require high-sensitivity data at resolu- 
tions below 1', that are pursued with current instrumen- 
tation and will eventually culminate with ALMA. 
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